arXiv:1504.01297vl [cond-mat.soft] 6 Apr 2015 


Rotation of a liquid crystal by the Casimir torque 


David A.T. Somers^’^ and Jeremy N. Munday^’^ 

^Department of Physics, University of Maryland, College Park, MD 20740, USA 
^Institute for Research in Electronics and Applied Physics, 

University of Maryland, College Park, MD 20740, USA 
^Department of Electrical and Computer Engineering, 

University of Maryland, College Park, MD 20740, USA 
(Dated: April 7, 2015) 

We present a calculation of the Casimir torque acting on a liquid crystal near a birefringent 
crystal. In this system, a liquid crystal bulk is uniformly aligned at one surface and is twisted at 
the other surface by a birefringent crystal, e.g. barium titanate. The liquid crystal is separated 
from the solid crystal by an isotropic, transparent material such as Si02. By varying the thickness 
of the deposited layer, we can observe the effect of retardation on the torque (which differentiates it 
from the close-range van der Waals torque). We find that a barium titanate slab would cause 5CB 
(4-cyano-4'-pentylbiphenyl) liquid crystal to rotate by 10° through its bulk when separated by 35 
nm of Si02. The optical technique for measuring this twist is also outlined. 


I. INTRODUCTION 

In 1948, Casimir [I] calculated that quantum fluctua¬ 
tions of electromagnetic fields cause attraction between 
two parallel metal plates at zero temperature. This can 
be interpreted as a physical manifestation of the zero- 
point energy predicted by quantum field theory [2]. As 
the plates are brought near each other, the allowed modes 
between them are limited. Because each electromagnetic 
mode has an energy of ?ia;/2 at T = 0 K, the total en¬ 
ergy of the configuration depends on the plate separation. 
This energy landscape results in an attractive force be¬ 
tween the plates. 

While Casimir’s description holds for ideal metal 
plates, a more general expression was derived by Lifshitz, 
Dzyaloshinskii, and Pitaevskii [3, 4]. This formulation is 
used for comparison with experiments and is valid for 
both real metals and dielectrics. In this theory, the force 
between two uncharged surfaces can be derived accord¬ 
ing to an analytical formula (often called the Lifshitz or 
Casimir-Lifshitz equation) that relates the zero-point en¬ 
ergy to the dielectric functions of the interacting surfaces 
and of the medium in which they are immersed. In 1972, 
Parsegian and Weiss [5] derived an expression for the non- 
retarded (where the speed of light is taken to be infinite) 
interaction energy between two dielectrically anisotropic 
plates immersed in a third anisotropic material. Barash 
[6] analyzed a similar problem including retardation ef¬ 
fects and found an equation for the Helmholtz free en¬ 
ergy per unit area. In the nonretarded limit, the results 
of Parsegian, Weiss, and Barash are in agreement. 

In the expressions derived by Parsegian, Weiss, and 
Barash, the Casimir energy between two parallel, bire¬ 
fringent plates depends on both their separation and their 
orientation. Between two parallel, positive birefringent 
slabs with in-plane optical anisotropy, the lowest energy 
state is when the two optical axes are aligned. This re¬ 
sults in a torque that rotates the plates to this configu¬ 
ration. 

While the Casimir-Lifshitz force is the subject of much 


discussion and has been verified in a number of experi¬ 
ments [7-12], there have been no published experimental 
attempts to measure the torque between anisotropic ma¬ 
terials. In addition to the early work of Parsegian, Weiss 
and Barash, recent theoretical work has been performed 
including a derivation of a more simplified equation of 
the torque between two plates in one dimension [13], 
the torque between two dielectric slabs with different di¬ 
rections of conductivity [14], and numerical calculations 
based on real materials [15, 16]. 

Most discussions of potential Casimir torque measure¬ 
ments involve either levitating microdisks [15, 16] or a 
torsion pendulum [17]. In the first method, microscale 
birefringent disks would be levitated over a birefringent 
substrate with Casimir-Lifshitz repulsion or an electro¬ 
static force. The Casimir torque would then cause the 
freely rotating disks to rotate so that their optical axes 
would align with that of the substrate. In the proposed 
torsion pendulum experiment, a macroscopic birefringent 
crystal would be attached to a quartz torsion pendu¬ 
lum and brought near another birefringent crystal. The 
Casimir torque would affect the period of natural oscilla¬ 
tions of the torsion pendulum, allowing for its detection. 

Here we propose a method that is in analogy to a static 
torsion pendulum with a thick liquid crystal layer as the 
twisted bulk. As the uniformly aligned liquid crystal is 
brought near a birefringent crystal, the Casimir torque 
aligns the liquid crystal molecules with the solid crystal’s 
optic axis, which in turn causes a twist through the bulk 
of the liquid crystal. A similar experiment was proposed 
by Smith and Ninham [18] in 1973 but, to our knowledge, 
was never carried out. Here we provide a calculation in¬ 
cluding retardation effects of the expected results of a 
similar geometry, as well as detailed experimental con¬ 
siderations. 

Instead of using the liquid crystal in the isotropic phase 
as a spacer layer (which is experimentally unfeasible), a 
thin layer of Si 02 separates the liquid crystal from the 
solid crystal. Varying the Si02 thickness is equivalent to 
changing the distance between parallel plates. If both the 
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FIG. 1. Experimental setup. The liquid crystal director is 
fixed at 45° at the top surface, but a Casimir torque at the 
bottom surface causes a linear twist throughout the bulk. In¬ 
cident light polarized at 45° is adiabatically twisted with the 
liquid crystal director to a final polarization state, 9f, which 
can be measured optically. Stronger Casimir torques cause a 
greater director twist. The inset shows a flat projection of the 
x-y plane. 


liquid crystal and the birefringent crystal have positive 
uniaxial anisotropy, then the Casimir torque causes the 
liquid crystal molecules to twist towards the extraordi¬ 
nary axis of the solid crystal. Because the liquid crystal 
is anchored at the glass interface, the director (i.e. the lo¬ 
cally preferred molecular orientation) is twisted through 
the bulk by the Casimir torque at the boundary at the 
opposite interface. This geometry is depicted in Fig. 1. 
This experimental design is similar to methods for mea¬ 
suring the azimuthal anchoring energy of liquid crystals 
on treated substrates [19]. The mechanism of liquid crys¬ 
tal alignment induced by a rubbed polymer layer itself 
is the subject of much study but is thought to include 
physical grooves, aligned polymer chains, and the van 
der Waals torque for surfaces in contact [20-23]. The 
last of these is equivalent to the short-ranged Casimir 
torque of our proposed study, though distinguishing the 
van der Waals effect from other alignment effects is diffi¬ 
cult. However, previous experiments have suggested the 
anisotropic van der Waals effect as the mechanism of liq¬ 
uid crystal alignment at a surface. Schadt and Schmitt 
[24] used linearly photopolymerized layers to align liquid 
crystals with a van der Waals interaction. Lu [25] also 
provided evidence that the van der Waals interaction is 
an important component of the liquid crystal alignment 
at treated polymer layers. However, these experiments do 
not isolate the van der Waals torque from other surface 
effects, because the liquid crystal is in contact with the 
substrate. Our proposed experiment would demonstrate 
this effect over a distance of tens of nanometers (and in 
doing so measure the long-range retardation effects of the 


FIG. 2. Drawing of the 5GB molecule, which is a nematic 
liquid crystal at room temperature. The extraordinary axis 
is along the molecular axis. Because 5CB has positive bire¬ 
fringence, Ue > Ho- 

Casimir torque). Finally, we can relate measured data to 
Casimir torques calculated from the dispersive properties 
of the materials. 

Smith and Ninham [18] considered the nonretarded 
case of this system and predicted measurable distortions 
of the liquid crystal director. Here we carry out the full 
retarded calculation of the Casimir torque by considering 
its effect on a thin boundary layer. That is, in comparison 
to the liquid crystal bulk with thickness t > 50 /im, most 
of the Casimir torque is felt by a thin layer of thickness 
5t < 50 nm. Also, because the total liquid crystal twist 
through the bulk is less than 45°, the liquid crystal in 
the region of 5t is nearly uniformly aligned. Therefore, 
we approximate the Casimir torque on this layer from 
the uniaxial crystal to be the same as that experienced 
by a uniformly aligned and semi-infinite liquid crystal 
slab. These approximations are treated with more detail 
in Sec. V. 


II. NEMATIC LIQUID CRYSTALS 

Most materials undergo a phase transition from solid 
to liquid at a melting point. Liquid crystals exhibit an 
intermediate phase in which there is some local ordering 
due to the molecular structure. Nematic liquid crystals 
are the simplest class. They can be thought of as long 
molecules (Fig. 2) that tend to align the long axis with 
their neighbors. The local direction of molecular orien¬ 
tation is written as a vector of unit length called the di¬ 
rector n. Because the molecules have different dispersive 
properties along the different axes, one can change the 
properties of the liquid crystal bulk by manipulating the 
director. We carry out calculations using the physical 
properties of 4-cyano-4'-pentylbiphenyl (5CB). In bulk, 
this is a positive uniaxial material so the dielectric func¬ 
tion is higher along the long molecular axis. 


III. CASIMIR INTERACTION BETWEEN TWO 
INFINITE SLABS 

To derive the interaction between two infinite slabs, 
Barash wrote the Helmhotz free energy between two bire¬ 
fringent materials at finite temperature as an infinite sum 
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TABLE I. Model parameters for dielectric functions of rele¬ 
vant materials. Oscillator data for BaTiOa, CaCOa, TiOa, 
and SiOa are from Ref. [31]. 




Cl 

oJi (eV) 

C2 

0J2 (eV) 

C3,4 

i^3,4 (eV) 

5CB 

T 

0.0374 

4.40 

0.1075 

5.91 

0.414 

9.19 


II 

0.0612 

4.40 

0.1025 

5.91 

0.460 

9.19 

BaTiOa 

T 

3595 

0.056 

4.128 

5.54 

- 

- 


II 

145.0 

0.138 

4.064 

5.90 

- 

- 

CaCOa 

T 

1.920 

0.138 

1.350 

13.4 

- 

- 


II 

1.960 

0.138 

1.377 

13.3 

- 

- 

Ti02 

T 

4.81 

5.069 

- 

- 

- 

- 


II 

5.62 

4.516 

- 

- 

- 

- 

Si02 


0.829 

0.057 

0.095 

0.099 

0.798 

1.098 

0.133 

13.39 


FIG. 3. Calculated torque between BaTiOs and uniform 5CB 
bulk separated by 10 nm of SiOa as a function of angle 0 be¬ 
tween the extraordinary axes of the birefringent materials. 
The line is a fit to a sin(20) dependence. The difference be¬ 
tween the fit and the full calculation is at most 2%, and is 
less than 0.1% for 0 — -k 


over Matsubara frequencies = {2'KkBT/h)n [6]: 
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where d represents the distance between the two slabs 
and 6 is the angular separation of their extraordinary 
axes. We have performed calculations for room tempera¬ 
ture based on experimental considerations; however, the 
zero-temperature calculation for the Casimir torque will 
provide the correct result to within <10% for the sepa¬ 
rations presented in this work (d < 50 nm). The zeros 
of D„(p,r) as a function of real frequency oj indicate al¬ 
lowed surface modes between the two infinite slabs. It 
can also be expressed in terms of the reflection matrices 
of the two materials, as in Refs. [26, 27]. In our regime of 
interest, the energy has sin^ 6 dependence to an excellent 
approximation, as in Fig. 3. The Casimir torque is then 


M{d, 9) 


dn 
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( 2 ) 


where a{d) is a negative for the materials considered here 
(which have positive birefringence). The dielectric func¬ 
tions of the two materials and of the intervening medium 
are evaluated at imaginary frequencies in Eq. (1). We 
use the Ninham-Parsegian oscillator model to describe 
the dispersion of the solid crystals [28]: 



d (nm) 


FIG. 4. Casimir torque per unit area between a half-slab of 
aligned 5CB and various birefringent crystals, when separated 
by a Si02 layer of thickness d with a relative angle of 45° 
between the extraordinary axes. 


with a three-oscillator model, so the dielectric function 
is: 


e5CB(*0 


l-b2 


3 

E 

i=i 


C, 


1-b 



( 4 ) 


For the birefringent materials, there are separate func¬ 
tions describing the ordinary and extraordinary axes. 
The model data used for our calculations are summa¬ 
rized in Table I. 


N 

*=(*0 = 1 +E I E ’ 

i=i ^ 

For the 5CB liquid crystal, we use the dispersive proper¬ 
ties at 298.2 K calculated by Kornilovitch [29] using data 
from Wu et al. [30]. There, the index of refraction is fit 


We calculate the Casimir torque between an infinite 
half-space of 5CB liquid crystal and an infinite half-space 
of several birefringent crystals using Eqs. (1) and (2) and 
the parameter data in Table I. The results for 9 = 45° 
are shown in Fig. 4. Note that torques per unit area on 
the order of 10 nN/m are found for separations of d = 10 
nm. 
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IV. TORQUE BALANCE METHOD FOR 
MEASURING LIQUID CRYSTAL ANCHORING 


The Casimir torque causing a director twist at one 
boundary competes with the restoring torque from the 
twisted liquid crystal. The latter is modeled using the 
Frank free energy density [32]. In our geometry, the di¬ 
rector is always aligned in the x-y plane so its orien¬ 
tation can be written in Cartesian coordinates as n = 
{cosd(z),sind(z),0}. There is no bend or splay of the 
liquid crystal, so only the twist term contributes to the 
distortion energy. The twist contribution is given by: 

Fd = ^{n-V xnf, (5) 


where k 22 = 3.6 pN is the twist elastic constant of the 
5CB[33]. Substituting our expression for n, the Frank 
free energy density is 


Fd = 


^22 



( 6 ) 


and the Frank free energy per unit area is: 

E.,...,. = (1) d,, (7) 

where t is the thickness of the liquid crystal layer. In our 
geometry (as in Fig. I), the extraordinary axis of the 
solid crystal is along the cc-axis, so the top boundary at 
z = d+t (where the liquid crystal is in contact with glass) 
is treated to induce alignment along 9{t + d) = 7r/4. 

A torque applied at z = d twists the director to 9{d) = 
Of. Using calculus of variations, the lowest energy conhg- 
uration is given by a linear twist, 0(z) = j + ^{z — d—t), 
where AO = -k jd, — 0 f . The elastic energy of the bulk per 
unit area is then: 


^elastic 


fez 
2 t 


( 8 ) 


If the director at z = d is twisted to d/, there is an energy 
penalty and associated restoring torque at that boundary. 
The restoring torque of the elastic bulk is given by: 


Alelastic 


^Aglastic _ /CZzAd 

dAO ~ i 


(9) 


A/eiastic is the torque that must be applied at the bound¬ 
ary z = d to twist the director to d(d) = Of. The torque 
applied at the boundary twists the director until the 
torque balance equation is satisfied: Meiastic+Afextemai = 
0. If the Casimir torque is in the approximate form 
Aleiastic ~ a{d) sin(2d), where 9 is the angle between the 
two extraordinary axes of the birefringent materials, then 
the torque is approximately Mcasimir = a(d)sin(2d/) = 
a(d) cos(2Ad) (which has the same form as the planar 
Rapini-Papoular approximation [34]), and the torque bal¬ 
ance equation yields: 


fez Ad 


-I- a(d) cos (2Ad) = 0. 


( 10 ) 


To predict the director twist for our proposed experi¬ 
ment, we calculate a(d) using Eqs. (1) and (2) and then 
numerically solve Eq. (10) to Hnd the twist caused by 
the Casimir torque. Ad. In the proposed experiment. Ad 
(or Of) will be measured to obtain a(d) via 

a(d) = ^^Adsec(2Ad). (11) 


V. THE BOUNDARY LAYER 
APPROXIMATION 

The liquid crystal can be treated as an anisotropic bulk 
material for the calculation because the twisting is small 
throughout the thickness that experiences the Casimir 
interaction. The spacing between the liquid crystal and 
solid birefringent crystal is on the order of ^10 nm, and 
the liquid crystal layer is about 100 /xm thick. Following 
Parsegian [35], the penetration depth of the Casimir in¬ 
teraction is on the order of the material separation. So, 
the most important region of the liquid crystal is the 10 
nm in contact with the SiOz (or conservatively, the 100 
nm). This 100 nm is the 0.1% of the liquid crystal near¬ 
est the birefringent crystal. Because the liquid crystal 
will be twisted a maximum of 7r/4 radians throughout 
the bulk, the liquid crystal director will vary by a max¬ 
imum of ~0.05° in the relevant region for the Casimir 
torque, which has no appreciable effect on the torque’s 
magnitude. 

This approximation can also be justified by consid¬ 
ering the reflection matrix of the liquid crystal stack. 
The Casimir energy of the system is a function of the 
reflection matrices of the two materials at the Matsub- 
ara frequencies. Our method assumes that, at the Mat- 
subara frequencies, the reflection matrix of the slowly 
twisted liquid crystal is nearly the same as that of an 
untwisted, bulk liquid crystal with the same alignment 
at the boundary. The first Matsubara frequency at room 
temperature, fe = 2TTkBT/h « 245 THz, corresponds to 
a wavelength of Ai « 1 pm. The higher frequences corre¬ 
spond to shorter wavelengths, so this first term has the 
longest penetration depth. We calculated the reflection 
matrices at this frequency between the twisted and un¬ 
twisted stacks using the Berreman 4x4 matrix method, 
and found them to be numerically identical to four sig¬ 
nificant Hgures [36]. The Casimir interaction energy is 
largely unaffected by the slow twist of the liquid crystal 
throughout the bulk. Hence, to several signihcant Hgures, 
the torque experienced by the liquid crystal layer is felt 
entirely at the nearest boundary and is only a function 
of the director orientation at that boundary. 


VI. PROPOSED EXPERIMENT 

Common methods for fabricating single liquid crystal 
cells have been previously reported in Refs. [37, 38] and 
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can be used for this experiment. A rubbed alignment 
layer of polyvinyl alcohol (PVA) can be used to cause 
the liquid crystal molecules to align along the rubbed 
direction at the surface. The birefringent crystal with a 
thin, isotropic Si 02 layer (with thickness d ^ 20 nm) is 
sandwiched with the PVA-treated glass with a spacing of 
t « 100 /i m (this value can be measured optically). The 
liquid crystal is then introduced into the cell via capillary 
action. The filling process may induce some alignment 
along the direction of liquid crystal flow; however, baking 
the sample above the LC clearing temperature (35° C for 
5CB) and allowing it to cool slowly will eliminate this 
effect. As the liquid crystal cools to room temperature, 
the director settles into the lowest energy state described 
in Sec. IV. The magnitude of the Casimir torque effect 
can then be measured by observing the twist of the liquid 
crystal director. 

The final director twist Of can be measured optically. 
This method is similar to a technique for measuring az¬ 
imuthal surface anchoring strengths of liquid crystals 
[39, 40]. When linearly polarized light is incident on 
an adiabatically twisted nematic liquid crystal stack (in 
which the pitch of the twist is much larger than the wave¬ 
length of light), the polarization state is rotated to follow 
the liquid crystal director. This is known as the adiabatic 
approximation for twisted nematics and is the principle 
behind twisted nematic liquid crystal displays [36]. In 
this experiment, white light polarized at 45° shines onto 
the stack as in Fig. 1. The Jones vector of this light 



, where I is the intensity. The light polar¬ 


ization follows the director of the twisted nematic and is 
incident on the transparent Si 02 layer with polarization 


df- 


Its Jones vector is now 



In a typical 


measurement of an anchoring force, the liquid crystal is 
sandwiched between two glass slides that do not interfere 
with the polarization state of the light, as in Ref. [19]. 
Then, the polarization state Of can be measured with a 
second polarizer. 

To calculate the expected results from an experiment, 
we consider a liquid crystal layer thickness of 100 /im and 
calculate the lowest energy state of the system using Eq. 
(10). Figure 5 shows the expected results for this case. 


The liquid crystal bulk is predicted to twist by over 35° 
when the stack is separated by 5 nm from BaTiOs, and 
a twist of several degrees is expected for separations of 
d ~ 50 nm (well into the Casimir regime). When near 
a birefringent material that has negative birefringence 
over a large frequency range (such as lithium niobate), 
the liquid crystal would twist towards the ordinary axis 
instead of the extraordinary axis. This would provide 
further confirmation that dispersion effects are causing 
the director to twist. 



d (nm) 

FIG. 5. Caculated twist of a 100 pm layer of 5CB caused by a 
Casimir torque induced by various birefringent crystals at dis¬ 
tance d from one end of the liquid crystal stack. The incident 
light is polarized at 45° to the ordinary axis at the top of the 
5CB stack, but the director is twisted by the Casimir torque, 
which causes the light polarization to rotate A9 towards the 
extraordinary axis at 0°. 


VII. CONCLUSIONS 


We have proposed an experiment for measuring a 
Casimir torque between a birefringent crystal and a liq¬ 
uid crystal separated by an isotropic spacer layer. We 
provide complete calculations of the expected results for 
several materials at a range of separations and include 
details for a proposed experiment. This experimental de¬ 
sign avoids many of the difficulties involved with a torsion 
pendulum or levitating microdisks. 
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